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' SUMMARY 

The variations of pressure distributions and coefficients of 
wave drag with fineness ratio, Mach number, and area ratio were 
calculated from linearized theory for a variety of open-nose 
bodies of revolution at zero angle of attack. Pressure distri- 
butions and values of wave-drag coefficients are given for sev- 
eral bodies with curved and straight lips and- for design and off- 
design conditions including the effects of additive drag. 

A comparison between the source method and the linearized 
method of characteristics showed that, for a given accuracy, com- 
putations by the source method were considerably faster; conse- 
quently, this method was used to compute the pressure distribu- 
tions for all bodies investigated. 


INTRODUCTION 

• Although methods for computing the supersonic flow over open- 
nose bodies of revolution have been known for some time (refer- 
ences 1 and 2), application of these methods has been limited to 
isolated cases. An investigation was therefore undertaken at the 
NACA Lewis laboratory to study systematically the variations of 
external pressure distributions and wave drags with fineness 
ratio, Mach number, and area ratio. The results presented herein 
are intended to serve as a guide in the design of external con- 
tours far supersonic nose inlets and to provide a basis for 
estimating the wave drag to be expected from such inlets . An 
evaluation of the relative merits of several methods of computa- 
tion is also included. 
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!Ehe following symbols are used in. this report: 
local area of cowling 

ratio of free- stream inlet area to maximum area 

free -stream Telocity of sound 

coefficient of wave drag, drag/q^AQ 

pressure coefficient, p-p 0 /q. 

maximum diameter of body 

functions of x/pR^ 

length of body 

fineness ratio 

free-stream Mach number, Tj/c 
static pressure 

dynamic pressure, pW 2 

radius of. body 
free-stream velocity 
axial-velocity increment 
radial-velocity increment 
total velocity 
cylindrical coordinates 


% 

P 


Mach angle, 
cotangent of 
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7 ratio of specific heats of air, 1.4 

€ Inclination of conical cowling with respect to free- 

stream direction 

c* transformed cowling angle, P tan e = tan €* 

» 

0 flow angle 

e c half -angle of cone 

0^ cowling positioning angle, angle between axLs of in l et 

* and straight line connecting vertex of cone with lip 

of cowling 

9 q shock angle 

p density 

Subscripts: 

0 free-stream conditions 

3 station of maximum diameter 

A,B,C,E characteristic grid points 
H nose of body 


BODIES AND VARIABLES INVESTIGATED 

The bodies investigated are shown in figure 1. Most of the 
computations were made for conical cowlings (fig.. 1(a)) to obtain 
the effect of free-stream Mach number Mg, area ratio Ag/Aj, and 
fineness ratio L/d. The ranges of the variables investigated for 
the conical cowlings were: Mach number, 1.5 to 2.5; area 

ratio, 0.400 to 0.800; and fineness ratio, 4 to 10. The effect 
of body contour on the pressure distribution and the wave drag 
for a Mach number of 2.0 was determined by investigating the flow 
over the curved cowlings (fig. 1(b)). Bodies for which the shock 
wave from the central body passes upstream of the cowling lip 
(fig. 1(c)) were considered in order to obtain the effect of addi- 
tive drag far a free-stream Mach number of 2.0. (The additive 
drag, as defined by Eerri and Hucci, is the integral of the pres- 
sure coefficient along the bounding streamline between the conical 
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shock and the cowling lip.) Pressure distributions were also cal- 
culated as a function of Mach number for a conical cowling with a 
central body. This central body was so located that the conical 
shock intersected the upstream edge of the cowling at a Mach num- 
ber of 2.0. 


COMPUTATIONAL METHODS 

Three of several methods available for calculating the pres- 
sure distributions over open-nose bodies of revolution are con- 
sidered: (l) the linearized source-distribution method (refer- 

ence l); (2) the rotational method of characteristics (reference 2); 
and (3) the method of linearized characteristics (suggested in 
references 3 and 4). A method applicable to conical cowlings has 
also been developed in reference 5. The method was not included 
in this Investigation, however, because in addition to the usual 
restriction to small slopes, the method requires that the cowling . 
be nearly cylindrical. 

The method of line arized characteristics, for which the equa- 
tions are derived in appendix A, will converge to the result 
given by the source method but the computing time for a specified 
accuracy is considerably greater (appendix B) j therefore the 
source method was used to compute the remaining pressure distri- 
butions. The rotational method of characteristics is the most 
accurate but the computations are more lengthy than those required 
for the linearized methods. The linearized theories give satis- 
factory results except for the cases where the Mach number or the 
body slopes are large. 

The wave-drag coefficient of each inlet was obtained by 
graphically integrating the pressure-coefficient distribution over 
the body surface. This coefficient was based on the free- stream 
inlet area and was defined by 
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where 



and A is the local area of the cowling. 


RESULTS AND DISCUSSION 

Conical cowlings . - The variations of pressure distribution 
in the axial direction for the conical cowlings are given in fig- 
ure 2 for selected values of the cowling angle and free- stream 
Mach number. In each case, the pressure coefficient has the 
Aokeret value at the lip, decreases rather rapidly, and at large 
distances downstream becomes asymptotic to the surface value pre- 
dicted for a cone of half -angle e . The variations of pressure 
distribution with fineness ratio and area ratio may be obtained 
by the correlation of the cowling angle with area ratio and fine- 
ness ratio (fig. 3). .The local pressures increase with increasing 
Mach number and decrease with increasing fineness ratio or area 
ratio. 

If a first-order approximation for the pressure coefficient 


°P 


P-PQ 

*0 



is used and if an appropriate change of variable in the axially 
symmetric, linearized potential equation is made, the quantity pCp 
becomes Independent of Mach number. By use of this approximation 
the results in figure 4 are obtained, where p 2 Cp is given as a 
function of the axial station x/pR^ for various values of the 
transformed cowling angle e’. 


Irr obtaining the curves presented in figure 4, the quan- 
tity p Cp was found to vary almost linearly with' e* for small 
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values ccf x/£%. For large? values of x/p% the variation 
of P 2 Cp with e • approached that of the cone solution, that 

is, (e‘) 2 log e This behavior suggested an attempt to fit the 
curves with an equation of the form 

t 

p 2 C p = f x e» - f 2 (e’) 2 logs e* (1) 

where e* is given in radians and fp and f« are functions 
of The values of fj_ and f g at x/gHu equal to zero 

can he immediately determined from the fact that the pressure coef- 
ficient at the lip is equal to the Acker et value, C p = 2e/p. 

Within the limitations cf linearized theory and with 
8 tan € = tan e’, the expression for the pressure coefficient 
at x/p% = 0 becomes =. 2e ' . Thus at the cowling lip, 

fp must equal 2 and fg must equal 0. When the remaining val- 
ues of f]_ and f 2 were obtained by matching the curves of 3 2 C p 
plotted against e’ at two points, equation (1) described the 
curves of P 2 C p plotted against e‘ accurately. As indicated 
by the preceding discussion, the term containing fp in equa- 
tion (1) is the predominating term near the lip and the term con- 
taining ±2 predominates at large values of x/pEjj. This fact 
is shown in figure 5 where the functions f p and fg are plotted 
for values of x/fS% from 0 to 10. It is expected that equa- 
tion (1) may be used over the range of e‘ for which the line- 
arized theory is valid. By use of equation (1), the general 
expression for the coefficient of wave drag of conical cowlings 
becomes 

/(% 

(fp-fgC* logg e*)(l+c*t) dt (2) 

where t = x/fffig. This expression was not used to compute the 
coefficients of wave drag because pressure distributions were 
separately calculated. If pressure distributions are not avail- 
able, however, the coefficient of wave drag can be obtained 
directly from equation (2). 

By integrating the pressure distributions of figure 2 over 
the body surfaces, the variations cf wave-drag coefficient with 
fineness ratio, Mach number, and area ratio are obtained 

(figs. 6 to 8, respectively). The coefficient of wave drag * 

decreases as fineness ratio, Mach number, or area ratio increases. 
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The "variation of the coefficient of wave drag with Mach number is 
relatively small compared with the variation with fineness ratio 
or area ratio. 

Although the wave drag decreases with an inc rease in fine- 
ness ratio, the friction drag increases. For large fineness ratios, 
the friction drag "becomes the predominant part of the total drag; 
hence to obtain an optimum cowling this fact must be taken into 
account. 

Additive drag and off -design conditions . - The coordinates »na 
contours of the inlets used to investigate the effects of additive 
drag are given in figure 9. The curves were arbitrarily chosen, 
except that the lip of each body is faired into a streamline of 
the conical field produced by the central body. The design param- 
eters for the bodies are 


Free-stream Mach number, Mq 2.0 

Area ratio, Aq/Aj 0.539 

Fineness ratio, L/D ...... ..... 4.0 

Cone half -angle, 6 C , deg 25 


The linearized pressure distributions for the three inlets 
having cowling positioning angles of 32°, 36°, and 40° are shown 
in figure 10. The linearized coefficients of wave drag for these 
bodies, the coefficients of additive drag obtained from linearized 
theory, and the coefficients of additive drag obtained from exact 
conical flow theory are given in the following table with the 
coefficients of additive drag and body drag for a cylindrical 
body (©2. = 30.5°) and a conical body (0^ - 42.6°): 


Cowling 

positioning 

angle, 

Coefficient of 
wave drag, Cp 
(body) 

Coefficient of additive drag 
CD. a 

Linearized cone 
theory 

Exact cone 
theory 

30.5 

0 


0.381 

32 

.016 

0.134 

.257 

36 

.029 

.046 

.088 

40 

.018 

.012 

.020 

42.6 

.014 

0 

0 


The linearized coefficients of wave drag for the bodies were com- 
puted by considering the bounding streamline between the conical 
shock and the cowling lip as part of the bodies; consequently, these 
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coefficients may "be somewhat in error because the slope at any 
point of the hounding streamline is large. The coefficient of 
total drag (not including a skin-friction coefficient) is the Stan 
of the coefficients of body and additive drag. Comparison of the 
coefficient of maximum total drag (0j = 30. 5°, cylindrical cowl- 
ing) with the coefficient of minimum total drag (07 *= 42.6°, 
conical cowling) shows that in this case the coefficient of total 
wave drag may be reduced by approximately 96 percent by eliminating 
the additive drag; therefore, as pointed out by I’erri and Eucci, 
the avoidance of this additive drag over a required range of free- 
stream Mach numbers is important in cowling design. 

Illustrated in figure 11 are the variations of pressure dis- 
tribution with Mach number for a cowling designed for a free- 
stream Mach number of 2.0. The area ratio of the conical cowling 
is 0.539, and the fineness ratio is 4.0. At the design Mach num- 
ber, the conical central body with a half -angle of 25° produces a 
shock that intersects the lip. For each off -design Mach number, 
the inlet flow is supersonic inasmuch as the normal shock is 
assumed to remain inside the inlet. The cowling in each case 
starts at the axial station x/Rjj of 0.159. The pressure dis- 
tribution upstream of this point for a free- stream Mach number 
of 1.8 results in an additive drag, because for any Mach number 
below 2.0 the conical shock is upstream of the cowling lip. As 
the Mach number increases, the pressure ratios on the cowling 
become higher. 

When the pressure distributions are integrated over the body 
to obtain the coefficients of wave drag and are plotted as a func- 
tion of free- stream Mach number (fig. 12), it is found that the 
coefficient of wave drag decreases rapidly from its value at a 
Mach number of 1.8 to its value at a Mach number of 2.0 and 
decreases slowly for Mach numbers higher than. 2.0. The coef- 
ficient of wave drag at a free-stream Mach number of 1.8 is 
approximately twice as great as that at a Mach number of 2.0; how- 
ever, again the slope at any point of the bounding streaml i ne is 
sufficiently large so that the drag coefficient obtained at a 
Mach number of 1.8 may be in error. For this reason, the accuracy 
of the dashed portion oif the curve in figure 12 between Mach num- 
bers of 1.8 and 2.0 is doubtful. 


Eonconical bodies without additive drag . - The coord in ates 
and the contours for the nonconi cal inlets are shown in figure 13. 
The design parameters are v 

Free-stream Mach number, Mq ............... 2.0 

Area ratio, Aq/a$ 0.539 

Fineness ratio, L/D • * 4.0 

Cone half -angle, 0 C , deg 25 
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The Initial lip angle for contours 1, 2, and 3 is the same as the 
flow-deflection angle behind the conical shock emanating from the 
central body. The results may again be somewhat in error because 
of the large slopes of the lips; however, the linearized theory 
should predict the drag trends of these bodies. 

Pressure distributions for the curved contours of figure 13 
are compared with the oonical -cowling pressure distribution in 
figure 14; the coefficients of wave drag for the four cowlings 
are presented in the following table: 


Contour 

Coefficient 
of wave 
drag 

°D 

1 

0.089 

2 

.042 

3 

.040 

4 

.014 


These values of the coefficient of wave' drag changed appreciably 
with the external contour of the cowling and the desirability of 
producing a rapid flow expansion near the lip is indicated. 


SUMMARY OF RESULTS 

The computations of the theoretical supersonic wave drag of 
axially symmetric open-nose bodies at zero angle of attack may be 
summarized as follows: 

1. For conical cowlings, the local pressures on a body 
increased with increasing Mach number and decreased with increas- 
ing fineness ratio and with area ratio. The coefficients of 
wave drag obtained from the pressure distributions decreased with 
increasing fineness ratio, Mach number, or area ratio. 

2. In order to determine the effect of additive drag, a 
series of cowlings having a common central body and the same mass 
flows were investigated at -a Mach number of 2.0. For these 
inlets, the least total drag was obtained with a conical cowl- 
ing the leading edge of which was located at the shock from the 
central body. As the leading edge was moved downstream from this 
location, the drag of the cowling decreased, but the increase in 
additive drag was sufficient to overbalan c e this decrease. 


0 
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3. For off -design, conditions, the total toys drag increased 
rapidly as the shock from the central "body passed ppstream of the 
cowling lip. "When this shock fell inside the lip, little change 
in wave-drag coefficient -with Mach number -was obtained. 

4. At a Mach number of 2.0, the coefficients of ware drag 
obtained for the nonconi cal bodies without additive drag varied 
from 0.089 to 0.040. The magnitude of the drag depended on the 
rate of flow expansion in the vicinity of the cowling lip. For 
an immediate flow exp ans ion, that is, a conical cowling, the coef- 
ficient of wave drag was 0.014. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, December 15, 1949. 


0 
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APPENDIX A 


33ERI7ATION OF EQUATIONS FOE LUCEAEIZED 
METHOD OF CHARACTERISTICS 

The linearized differential equation for supersonic flow Is 
of Hyperbolic form and therefore has real characteristics. The 
use of these characteristics for the numerical computation of 
supersonic flow fields when the perturbation velocities are small 
has been suggested by several authors. (See, for example, refer- 
ences 3 and 4.) Inasmuch as these characteristic equations are 
very simple to use for calculating axially symmetric flow, a com- 
parison of this method with the source method was considered 
desirable. The required difference equations Eire derived herein. 

For axially symmetric flow, the linearized equation for the 
perturbation velocities u and v is 

(Ai) 

dx Sr r 

The characteristics of this equation are given by (reference 4) 

dr 1 

dx p 

so that the total derivatives cf the perturbation-velocity com- 
ponents became 


du = du 1 du 
dx dx p Sr 

(A2) 

dv = civ ± 1 dv 
dx dx p dr 

(A3) 


where the plus sign is used for family I and the minus sign is 
used for family II (fig. 15(a)). By solving equations (A2) 
and (A3) for du/dx and dv/ar and by making use of the irro- 
tationality condition, equation (Al) for family I becomes 
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dU -^ T - -kz dz = 0 (A4) 

p r 

’When the grid shown In figure 15(a) is used, equation (A4) may he 
expressed in difference form as follows: 



Similarly, for family H, the difference equation is 

u c+j = li B+y + ^ ^ = f b 

or 

Uq = pT + (^6) 

By solving equations (A5) and (A6) simultaneously, u^ and v c 
are found to he 

n 0 - 1 <W (A7) 

T c - | » B - V (A8) 

The flow angle 0 and the total velocity W are introduced by 
letting 


u +TJ = W cos 0 


v = W sin 0 


(A9) 

(A10) 
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From equations (A7) to (A10), the expressions for 6 C and W c 
are found to he 


vhere 


e = ten-1 

° V°A 


v„ 


2 sin. 6q 


G A " V A 008 e A 


, tane A ^ 

e V P r A> 


Gg = W B cos e B 


1 + ^A 


P 


Pr- 


B 


(All) 

(A12) 


By considering figure 15(h), the oppressions for the coordinates 
and r c may he derived as 


*C =! (*B+*A> + | (^B~ r A) 


r c c r c - 




P = r A + 


^C-^A 

P 


(A13) 

(A14) 


lor the velocity components at a body boundary point 
(fig. 15(c)), the difference equation is 

+ T ‘ + T + p (?)o (A15) 

Ina smuch as the flow angle at the bo undar y is known, one of the 
unknown velocity components in equation (A15) can he eliminated 
with the aid of the boundary c ondi tion 

CA1B) 
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By making the substitutions given in equations (A9) and (A10), 
% is found to be 




cos 6q 
E ~ "C cos ©g 


V. 


tan Q r 


1 + 


1 + 


*E mX G \ 

) 


1 + 


tan 


(A17) 


With the velocity known, the static pressure at a boundary point 
is determined from 



(A18) 


The coordinates Xg^. and Bg for the boundary points are 




r C - r A ^ ~e~ + X A ^ e A 

I +tone A 


Bg - r c - 


P 


(A19) 


(A20) 


The linearized method of characteristics differs from the 
exact method (reference 2) in that all characteristics are straight 
and parallel so that the characteristic grid can be constructed 
in advance of any computations. All the coordinates of the grid 
points are therefore known end need not be computed if a graphical 
representation is employed. The equations for the velocity com- 
ponents (equations (A12) and (A17 ) ) are, of course, also consider- 
ably simpler than those obtained for the exact method of 
characteristics . 



NACA TN 2115 


15 


APEENDH B 


COMPARISON OF ACCURACY AND COMPUTING TIME 

A comparison of .the pressure distributions computed by the 
linearized method of characteristics and by the source method is 
shown in figure 16, One curve was obtained with a spacing 
between bou ndar y points of Ax/% = 0,5 and a second with 
Ax/% = 0,125, The third curve was obtained by the source method 
with 11 boundary points. It was found for the source-distribution 
method that any number of boundary points above 10 did not appre- 
ciably alter the pressure-distribution curve. The calculations 
by the method of linearized characteristics, as presented herein, 
require a larger number of boundary points than the source method. 
The asymptotic solution (corresponding to zero grid spacing), 
however, should be the same for the two methods inasmuch as both 
methods are derived from the same differential equation. 

A comparison of the source method with the method of rota- 
tional characteristics is presented in figure 17 and shows a 
relatively large error in the vicinity of the cowling lip. The 
wave-drag coefficients obtained by integrating the pressure dis- 
tributions of figure 17 to the point x/% = 2 are: exact- 

characteristics method, O.lOOj and source method, 0.064. The mag- 
nitude of the error is approximately 36 percent and should repre- 
sent a near maximum. Inasmuch as the initial flow deflection at 
the lip Is large (15.5°), this error is to be expected because 
linearized theory Is inaccurate for large angles. Some idea of 
the magnitude of the error due to the limitations of linearized 
theory (in regard to flow-deflection angles) can be obtained 
from figure 18, where the pressure . coefficients obtained from 
linearized theory and from shook theory are plotted as functions 
of the lip or initial flow-deflection angle. 

In making a calculation by the source method for a body 
having 15 boundary points (average number of points used in the 
calculations), the computing time, including the time required 
for checking the computations, was approximately 17 man-hours. 

If the number of boundary points was increased to 20, the time 
necessary for computing and checking would be 30 man-hours. The 
time to make the same calculation (curve with Ax/% = 0.5, 
fig. 16) by using the linearized method of characteristics was 
approximately 160 man-hours. This time, however, includes com- 
puting the coordinates of the grid points and could be shortened 
if these were read from the grid network. The same calculation 
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made "by the method of characteristics for rotational flow required 
approximately 240 man-hours. The over-all time for a calculation 
by any one of the preceding methods depends, of course, on the 
experience of the computers should not be considered a min i m u m 
value. 
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(a) Conical cowlings. 



(t>) Curved cowlings. 



(c) Cowlings with additive drag. 
Figure 1. - Bodies investigated. 



tic-pressure ratio 







Statio-presaure ratio, p/p 





Axial station* x/R 

xi 

(o) Free-stream Mach number, 2.5. 

Figure 2. - Concluded. Pressure distributions in axial direction for conical cowlings 
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Fineness ratio, L/D 


Figure 5. - Relations among cowling angle, fineness ratio, and 

area ratio® 
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Figure 4. - Pressure distributions obtained from first-order 
equation for pressure coefficient. P tan e = tan e 1 . 
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Figure 5. - Variation of functions of x/(3R^. 



(a) Free-stream Maoh 
number, 1.5, 

Figure 6. - Variation 


Fineness ratio, L/b 

(b) Free-stream Mach (c) Free-str 

number, 2.0. number, 

of coefficient of wave drag with fineness ratio. 


co m 

















Coeffi 






Contour 


Center line 


V R 0 


R 2 //R 0 


Rj/Rq 1.026 1.048 1.060 1.066 1.076 1.083 1^086 1.090 1.096 1.106 1.123 1.202 1.280 1.300 


0.956 

1.176 

1.260 

1.310 

1.190 

1.210 

1.090 

1.096 


1.360 1.360 1.360 1,560 1.560 1.560 


1.10011.13011.160 1.170 1.190 1.210 1,250 1.240 1.280 1.310 1.358 1.360 


Figure 9. - Cowling inletB need to Investigate effect of additive drag. Free-atreaia Hack number 
2.0$ area ratio, 0.639j fineness ratio, 4.0$ cone half-angle, 25°. 
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0 2 4 6 8 10 12 

Axial station, x/Rq' 


Figure 10. - Pressure' distribution along bounding streamline and* 
over body contour for inlets with additive drag. Free-stream 
Mach number, 2.0; area ratio, 0.539; fineness ratio, 4.0. 
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Figure 11. - Variation of pressure distribution with Mach number 
for a cowling designed for free-stream Mach number of 2.0. 

Area ratio, 0.539; fineness ratio, 4.0; cowling angle, 1.9°. 



1.8 2.2 2.6 3.0 


Free-stream Mach number. Mg 

Figure 12 c - Coefficient of wave drag as function of Mach number 
for nose inlet designed for free-stream Mach number of 2.0. 
Area ratio, 0.5395 fineness ratio, 4.0. 
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1.36 

1.36 

1.56 
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i 
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Figure 13. - Coordinates and contours for nonconical Inlets. Free-stream Mach number, 8.0$ 

area ratio, 0.639$ fineness ratio, 4.0. 
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Figure 14. - Pressure distributions over several body contours. 
Free-stream Mach number, 2.0; area ratio, 0.539; fineness ratio 

4„0. 
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Figure 16. - Comparison of pressure distributions obtained with, 
linearized characteristics and with source method. Free-stream 
Mach number, 1.5; area ratio, 0.393; fineness ratio, 4.0. 
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Figure 18. - Lip-pressure coefficient as function of lip angle. 
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